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ABSTRACT 


The  dynamic  processes  and  stability  of  plasma  discharges 
in  the  intermediate  range  between  the  spark  and  the  arc  type  of 
discharge  in  liquid  dielectrics  are  examined.  The  object  of  this 
investigation  is  the  industrial  application  of  this  type  of  dis¬ 
charge,  in  particular  electro-discharge  machining.  Available  in¬ 
formation  about  this  topic  is  assembled  and  discussed.  In  part¬ 
icular  the  low  frequency  liquid  oscillations,  their  damping  and 
thermal  instabilities  are  investigated. 

The  results  of  experiments  performed  to  investigate  the 
properties  of  discharges  in  water  and  other  liquid  dielectrics 
are  presented.  A  discussion  of  experimental  and  theoretical  re¬ 
sults  with  respect  to  practical  application  is  given. 
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CHAPTER  1 


INTRODUCTION 


1.1  Industrial  applications  of  plasma  discharges  in  liquids 

Plasma  discharges  have  been  employed  in  many  branches  of 
industry  for  many  years.  Probably  the  best  known  is  the  high  pressure 
arc  burning  in  a  gas.  The  high  temperature  of  its  plasma  column  lends 
itself  to  welding,  lighting  and  in  metallurgical  processes. 

Other  industrial  applications  for  the  electrical  discharge 
have  been  found  in  the  past  30  years.  In  particular,  discharges  in 
liquids  which  are  convenient  for  many  technological  operations  have 
been  used.  These  can  be  divided  in  three  main  groups: 

i)  electro-discharge  machining  (EDM) 

ii)  electro-hydraulical  effect  (EHE)  -  production  of  powerful  shock 
waves  whose  applications  include  the  forming  of  metal  plates  and  the 
disintegration  of  minerals. 

iii)  plasma  torches  -  used  for  metal  cutting,  surfacing,  chemical 
synthesis,  etc. 

At  present,  the  most  widely  used  technology  of  the  above 
three  'groups  is  EDM.  The  principle  of  this  method  consists  of  an 
errosion  effect  caused  by  the  spark  (or  not  purely  spark)  discharge 
on  the  electrode  surfaces  in  the  liquid.  The  discharge  burns  between 
two  electrodes  in  a  gap  about  0.05  to  0.2  mm  wide  normally  filled  with 
a  liquid  dielectric  such  as  water  or  oil.  Values  for  the  discharge 
duration  range  from  1  to  1000  psecs.  Many  changes  in  the  liquid  and 


. 
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on  the  electrode  surfaces  may  be  observed  after  the  discharge 
extinguishes.  The  most  important  change  is  the  creation  of  a  crater 
on  the  electrode  surface.  Under  certain  conditions  the  volume  ratio 
of  the  two  electrode  craters  may  extend  from  1:1  to  1:1000.  Through 
the  use  of  a  large  number  of  discharges  a  machining  process  is  obtained 
where  the  electrode  with  the  smaller  craters  is  a  tool  and  the  electrode 
with  the  larger  craters  is  the  workpiece. 

1.2  Spark  and  arc  discharges 

The  basic  difference  between  the  spark  and  the  arc  discharge 
lies  in  the  dynamics  of  the  two  phenomena.  The  arc  is  a  steady  plasma 
discharge  with  no  considerable  instabilities.  A  spark  on  the  other  hand 
is  a  very  dynamic  discharge;,  changing  its  geometry,  temperature, 
pressure  and  composition  on  a  very  short  term  scale,  typically  several 
hundred  nanoseconds.  A  spark  is  often  the  initial  state  of  the  arc 
discharge.  A  very  rapid  expansion  of  the  plasma  channel  is  a  typical 
characteristic  of  a  spark. 

The  velocity  of  this  expansion  is  greatest  in  the  initial 
state  and  reaches  values  of  several  hundred  meters  per  second.  The 
velocity  of  the  initial  expansion  of  the  plasma  channel  is  therefore 
supersonic  and  is  accompanied  by  one  or  more  shock  waves.  While  the 
overpressure  inside  the  arc  column  is  very  small  (few  centimeteis  of 
the  water  column) ,  the  spark  overpressure  is  very  high  and  can  reach 
several  ten  thousand  atm.  of  pressure  in  the  liquid  environment.  The 
high  overpressure  is  obviously  connected  with  different  temperatures 
inside  the  discharge  channel.  The  temperature  inside  the  spark  dis¬ 
charge  channel  is  considerably  higher  than  that  inside  the  arc  discharge 
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assuming  that  both  discharges  occur  under  equal  conditions. 

The  original  EDM  machines  employed  only  a  spark  discharge. 
The  spark  discharge  was  driven  from  a  so  called  relaxation  generator 
shown  in  Figure  (1.1). 


Fig.  (1.1)  The  working  principle  of  the  relaxation  EDM  generator. 

The  capacitor  C  is  charged  from  the  power  source  through 
resistor  R  and  is  discharged  into  the  spark  gap  after  the  voltage  on 
the  capacitor  is  sufficiently  high  for  liquid  breakdown  to  occur. 

The  low  inductive  connection  between  the  capacitor  C  and  the  spark 
gap  allows  a  very  short  high  amplitude  current  discharge. 

The  modern  generators  (so  called  independent  generators) 


with  more  flexible  parameters  are  able  to  provide  impulses  of  the  type 
shown  in  Figure  (1.2) 


- - t 

Fig.  (1.2)  The  working  principle  of  the  independent  EDM  generator. 
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By  using  a  triggered  switching  element,  current  impulses 
(usually  rectangularly  shaped)  with  an  amplitude  less  than  that  de¬ 
picted  in  fig.  (1.1)  may  be  obtained.  The  length  of  such  impulses 
may  reach  several  hundreds  psecs. 

The  independent  generators  do  not  work  with  purely  spark 
discharges.  The  supersonic  expansion  lasts  only  a  very  short  time 
at  the  very  beginning  of  the  discharge.  During  the  remainder  of  the 
discharge,  the  expansion  velocity  is  below  supersonic  and  the  high 
overpressure  typical  of  a  spark  no  longer  exists.  However,  the 
depression  and  sometimes  also  the  channel  contraction  not  typical  of 
an  arc-like  steady  discharge  may  still  exist. 

Therefore  this  kind  of  discharge  may  be  divided  in  three 
phases.  The  first  phase  is  the  spark-like  behaviour,  the  second  an 
intermediate  subsonic  expansion  and  the  third  an  arc— like  burning 
under  changing  outer  conditions. 

For  EDM,  the  first  and  third  phases  of  the  discharge  are 
of  most  importance.  The  first  phase  is  of  importance  because  of  the 
high  temperature  inside  the  plasma  channel  and  the  strong  shock  waves 
propagating  from  the  discharge.  The  third  phase  is  important  because 
the  discharge  is  in  this  phase  for  an  essential  part  of  the  discharge 
time. 

The  life  time  of  the  second  phase  is  relatively  short,  making 
this  portion  of  the  discharge  insignificant  in  the  study  of  EDM. 

Thus  the  investigation  of  the  independent  EDM  generator  dis¬ 
charges  may  be  separated  into  the  study  of  the  spark— like  and  the  arc 


like  phases. 
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The  first  spark-like  phase  in  which  a  relatively  small  amount 
of  energy  is  released  is  a  rather  complicated  problem  and  will  not  be 
discussed  in  this  thesis.  The  influence  of  the  first  and  second  phases 
on  the  third  will  be  discussed  only. 

1.3  Importance  of  discharge  instabilities  in  EDM 

Steady  discharges  with  the  released  energy  preserved  by  using 
the  minimum  possible  current  amplitude  are  suited  to  many  EDM  operations. 
However,  the  lowering  of  the  current  amplitude  is  limited  by  the  dis¬ 
charge  instabilities  which  are  considerably  enhanced  at  small  current 
amplitudes.  Any  changes  in  the  shape,  position  and  magnitude  of  the 
discharge  channel  is  for  technological  reasons  undesirable  because 
according  to  experimental  data  these  are  connected  with  the  tool 
electrode  damage.  This  damage  is  probably  caused  by  shock  or  other 
mechanical  waves  which  accompany  the  plasma  column  instabilities,  hit¬ 
ting  the  overheated  and  melted  electrode  surface.  This  results  in 
pieces  being  torn  out  of  the  electrode  material. 

It  is  known  that  under  certain  outer  conditions  it  is  possible 
to  operate  with  very  small  current  amplitudes  without  any  sensible  in¬ 
stabilities.  On  the  other  hand  with  some  different  technological  opera¬ 
tions,  instabilities  are  observed  with  much  higher  current  amplitudes. 
Thus  the  instabilities  may  be  influenced  by  outer  conditions.  The  in¬ 
stabilities  are  manifested  in  the  outer  electrical  circuitry  as  changes 
in  the  discharge  current. 

If  the  instabilities  are  sufficiently  intense  and  the  current 
amplitude  sufficiently  small,  discharge  interruption  can  occur.  From 
the  technological  point  of  view  this  is  the  worst  case,  because  each 
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interruption  implies  a  new  start  of  discharge  with  the  additional 
undesirable  spark-like  discharge  phase. 

A  complete  theoretical  description  of  the  type  of  discharge 
discussed  above  has  not  been  done  yet.  The  theory  describing  a  pure 
spark  discharge  is  satisfactorily  developed  and  is  successfully  applied 
in  the  relaxation  generator  EDM  theory.  The  theory  of  the  stable  high 
pressure  arc  discharge  is  also  well  developed.  The  results  of  the 
theory  on  pure  spark  discharges  and  the  theory  on  high  pressure  arc 
discharges  may  be  used  as  a  basis  for  developing  a  more  detailed 
theoretical  and  experimental  investigation  of  the  independent  EDM  dis¬ 
charges.  A  complete  study  of  the  discharge  instabilities  is  a  rather 
complex  problem.  An  instability  may  develop  for  several  reasons,  in 
addition  the  coexistence  of  several  kinds  of  instabilities  can  not  be 
excluded. 

In  addition  many  instabilities  as  described  in  publications, 

1  2 

still  await  a  theoretical  description,  eg.  Hoyaux  ,  Kocian  . 

The  important  unknowns  in  the  independent  generator  EDM 
method  are  the  magnitude  of  the  pressure  and  temperature  in  the  dis¬ 
charge  space  during  and  after  the  discharge  and  the  origin  and  state 
of  the  discharge  channel  instabilities.  These  problems  can  not  be 
exactly  theoretically  solved.  An  exact  theoretical  solution  for  these 
unknowns  is  not  possible  due  to  the  complexity  of  the  phenomena  involved. 
Some  simplification  in  the  theoretical  description  of  the  discharge  is 
necessary  if  an  approximate  solution  is  to  be  obtained.  An  experimental 
verification  of  the  simplifying  assumptions  must  then  be  carefully  made. 
Similarly  all  experimental  data  should  be  properly  theoretically  explained. 
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In  this  thesis,  contemporary  knowledge  about  important 
phenomena  occurring  in  plasma  discharges  is  collected.  Phenomena 
such  as  the  oscillation  of  the  gas  bubble  and  the  thermal  instability 
are  further  investigated.  The  results  of  an  experiment  performed  to 
investigate  the  discharge  instability  under  various  conditions  in  a 
liquid  dielectric  will  be  discussed  and  compared  with  theoretical 
results.  The  main  objective  of  the  thesis  is  to  establish  the  most 
important  phenomena  in  discharges  occurring  in  liquids. 


' 
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CHAPTER  2 


DYNAMICS  OF  THE  STABLE  DISCHARGE  IN  A  LIQUID 


2.1  Plasma  channel  created  by  a  spark  discharge  in  gases  and  liquids 

The  dynamics  of  a  spark  discharge  in  gas  or  liquid  have  been 

3  4  5 

discussed  by  many  authors  (Skvortsov  et  al  ,  Abramson  et  al  ’  , 

Drabkina^,  Gegechkori^,  Braginskii^,  Andreev  et  al^’^,  Egorova  et  aldd)  . 

All  authors  assume  a  cylindrical  shape  for  the  radially  expanding  plasma 

column  which  is  pushing  a  shock  wave  into  the  undisturbed  environment. 

The  expansion  of  a  spark  channel  in  a  gaseous  environment  was  first 

4-  6 

described  by  Drabkina  et  al  5  in  1945  and  1951. 

3 

According  to  Drabkina  the  discharge  energy  is  converted  into 

the  creation  of  a  pressure  disturbance  with  a  shock  on  the  edge  of  the 

3 

disturbance.  By  using  the  hydrodynamic  equation  Drabkina  was  able  to 
solve  for  the  distribution  of  the  pressure,  density  and  temperature 
behind  the  front  edge  of  the  expanding  disturbance.  The  following 
equations  describe  the  motion  of  the  radius  of  the  channel  and  the 


shock  wave 


rg(C) 


!/ 4 

po 


t  1/2 

f  E(t)  1/2 

<  t  >  dt 


o 


(2.1) 


PN^  E ( t )  ^ 
rk(t)  -  L(— ) 

o 


Em  1/2 

( i  }  dt 


o 


N 


(2.2) 


where  r  (t)  -  channel  radius,  r  (t)  -  shock  wave  radius,  t  -  time, 
s  k 

E(t)  -  released  energy  up  to  time  t,  £  -  length  of  the  spark,  - 
gas  density  at  760  torr,  p'  -  gas  density  in  the  observed  space, 
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a  ,  L,  p,  M,  N  -  constants  depending  on  the  gas  properties. 

The  Drabkina  theory  has  been  discussed  and  experimentally 

3  5  7 

verified  by  other  authors  .  In  1958  a  more  accurate  theory  was 

g 

developed  by  Braginskii  which  was  later  partially  improved  and 
completed.  According  to  many  authors,  the  Braginskii  theory  is  in 
very  good  agreement  with  present  experimental  data. 

In  the  Braginskii  theory  the  plasma  channel  is  assumed  to 
be  a  homogeneous  cylinder.  The  energy  released  in  the  spark  is 
accounted  for  by  an  increase  in  the  internal  energy  (dissociation, 
ionization,  increased  temperature)  and  by  the  expansion  of  the  gas 
in  the  channel. 

If  the  effects  of  radiation,  skin-effect  and  magnetic 
forces  are  neglected  and  the  electrical  conductivity  assumed  approxi¬ 
mately  constant  throughout  the  plasma  column  diameter,  the  quasi-self- 
similar  solution  for  the  cylindrical  shock  wave  from  the  energy  balance 
equation  gives  the  following  expression  for  the  channel  radius: 


rk(o  =■ 


2  _ 

TT  Pq50 


1/6 


J2/3dt 


1/2 


(2.3) 


o 


where  J  -  electric  current,  a  -  electric  conductivity,  p  -  constant 
depending  on  the  gas  properties. 

The  Braginskii  theory  has  been  successfully  used  by  many 
authors9, 10 ’ 11 .  The  best  discussion  and  comparison  between  Drabkina' s 

experimental  data  and  the  Braginskii  theory  appears  in  the  work  of 

1 9  13  12 

Koppitz  and  Tholl  .  Koppitz  has  improved  Braginskii  s  expression 
for  the  latter  stages  of  the  discharge  where  the  initial  assumptions 


are  no  longer  fulfilled. 
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For  the  description  of  a  spark  discharge  in  a  liquid  both  the 

g  5 

Braginskii  and  Drabkina  theories  may  be  used.  The  difference  between 
the  case  of  a  discharge  in  a  gaseous  medium  and  that  of  a  discharge  in 
a  liquid  is  that  in  the  case  of  a  liquid  the  surrounding  medium  is  less 
compressible  and  is  more  dense.  As  a  result  more  energy  is  converted 
into  the  creation  of  the  shock  wave  and  into  the  kinetic  energy  of  the 
moving  liquid.  Naturally  the  pressure  inside  such  a  discharge  is 
considerably  higher  than  in  a  similar  discharge  in  a  gaseous  environ¬ 
ment. 

The  cylindrical  model  will  also  be  used  to  describe  the  spark 
discharge  in  a  liquid  (see  fig.  2.1a).  The  assumed  density  and  temp¬ 
erature  profiles  in  the  gas  bubble  are  depicted  in  Figure  (2.1b). 


Fig. (2.1)  The  cylindrical  discharge  model  with  temperature 
T  and  density  p  distribution 

2.2  The  steady  plasma  column  of  a  discharge  in  a  liquid 

If  liquid  vaporization  is  neglected,  the  stable  arc  dis¬ 
charge  in  a  liquid  environment  is  similar  to  a  discharge  in  a  solid 
tube  with  a  variable  diameter.  Such  a  discharge  is  described  in 
many  publications,  the  most  important  of  which  are  by  Hoyaux\ 
Finklenburg"^ *  ^  Maecker^^  ,  Ecker  and  Zoller^. 

The  high  pressure,  high  current  and  the  thermally  in¬ 
homogeneous  plasma  column  is  described  in  detail  by  Ecker  and  Zoller^. 
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Their  theory  does  not  require  the  Shotky  assumption  of  a  constant 
temperature  across  the  plasma  column.  Their  results  provide  the 
basis  for  the  thermal  instability  described  by  Ecker,  Kroll,  Zoller 
The  theory  is  based  on  a  model  consisting  of  a  infinitely  long 
cylindrical  plasma  column  with  insulated  walls  and  a  radius  R .  The 
plasma  is  taken  as  being  weakly  ionized  across  the  column.  The 
result  of  the  analysis  is  the  following  system  of  two  coupled  differ¬ 
ential  equations  which  include  all  temperature  effects: 

~  -57  (rD[T]  +  i  ^  ( rnD [ T ]  In  T)  +  a[T]n  =  0  (2.4) 


r  df  (rA[T)  +  elj2  V11"  =  ° 


(2.5) 


with  the  boundary  conditions 

dn 


r  =  0 


r  =  R 


dr 


=  0 


dT 


=  0 


n  =  0 


dr 
T  =  T 


(2.6) 


wall 


where  D  -  ambipolar  diffusion  coefficient,  r  -  radius,  X  -  heat 
conductivity  coefficient,  T  -  neutral  gas  temperature,  a  -  particle 
production  coefficient,  n  -  particle  density,  e  -  particle  charge, 

E  -  longitudinal  electric  field,  y^  -  electron  mobility. 

These  equations  reveal  a  similarity  law.  They  are  not 
analytically  solvable  and  therefore  a  numerical  solution  is  necessary. 
The  similarity  law  is:  In  plasma  columns  with  identical  values  of 
Rp  and  E/p  the  radial  temperature  distributions  are  identical  and  the 
density  distributions  similar.  The  density  is  directly  and  the  current 


inversely  proportional  to  the  pressure. 


A  typical  change  of  the  radial  neutral  gas  temperature  and 


' 


■ 
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Fig.  (2.2a)  Typical  change  of  the  radial  gas 
temperature  distribution  (after  Ecker 
and  Zoller-*-^) 


Fig.  (2.2b)  Dependence  of  the  axial  gas  temperature 
T„  on  pressure  P,  current  I  and  radius  R 
(after  Ecker  and  Zoller  °) 
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the  dependence  of  the  axial  neutral  gas  temperature  for  a  helium 
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column  as  calculated  by  Ecker  and  Zoller  are  shown  in  Fig.  2.2a 
and  2.2b. 

2.3  Oscillation  of  the  gas  bubble 

For  the  further  investigation  of  the  processes  in  the  spark 
gap,  a  description  of  the  behaviour  of  a  simple  gas  bubble  in  the  gap 
area  is  useful.  In  particular,  when  the  spark  discharge  extinguishes 
the  remaining  gasses  build  up  just  such  a  bubble.  The  behavious  of 
the  bubble  is  determined  by  the  conditions  which  prevailed  at  the  end 
of  the  discharge.  In  other  words,  the  pressure,  temperature,  etc. 
which  prevailed  at  the  end  of  discharge  can  be  deduced  if  a  suitable 
record  and  description  of  the  gas  bubble  is  available. 

The  shape  of  the  bubble  depends  on  the  shape  of  the  electrodes. 
A  cylindrical  shape  for  the  bubble  will  be  assumed  for  the  case  of  two 
parallel  electrodes  and  a  spherically  shaped  bubble  for  the  case  of 
two  tip  electrodes.  The  first  approximation  to  the  solution  for  the 
oscillation  of  a  gas  bubble  will  be  obtained  by  considering  the  ideal 
case  with  no  friction,  no  heat  exchange  between  bubble  and  environment 
and  a  perfect  gas  inside  the  bubble. 


Fig.  (2.3)  The  oscillating  cylindrical  bubble 
The  equation  of  motion  for  a  cylindrical  bubble  in  Fig.  (2.3) 
can  be  derived  from  Bernoulli's  equation  as  follows: 
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. 
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where  c|> 
function 


for  r  = 

but 


where 

then  p  = 


but 


3<b  .  1  ’  2  .  P  r  /  \  '  '  R  *  3<j>  A  - 

3t  2  p  K  ’  *  r  »  3r  »  Y 

-  velocity  potencial,  p  -  liquid  density,  f(t)  -  arbitrary 

of  time,  t  -  time. 


P  ,  ,  3  (RR)  ,  1  ^2  R  c  f  i-\ 

-  +  la  r  -j-j—  +jR~2=  f(t) 


for  r  ->  03  f  (t)  =  0 


P  *2  i  .  2  R 

-  =  -[In  r  (RR  +  R;  +  ~  R  —  ] 
p  2  Z 

K  r 


l  we  get 


P  »»  •  2  1*2 

~  =  -[In  R  (RR  +  R  )  +  -  R  ] 
P  4 


A  j>2\  _  d  iR  it 
R  ’  dR  (2  R  )  ~  R  dt  dR 


t  =  +i2)  lnR  +  ii2' 


For  an  adiabatic  compression  and  for  r  ->  00 ;  p  ->  p*  we  get 


Pi-=(-r->2y 
'‘ID  / 


p  +  p*  R^ 


=  p*  +  pCRj^) 

R-,  2y 

pRr-) 


R^  -  minimal  radius  or  initial  radius 
y  -  specific  heat  ratio 


(R  I  Ir  +  r2)  ln  R  + 1  r2  -  f  -  r  (^)2Y 


.2  dR  .  -2, 


.2  2 


R. 


2y 


(R**  ~  +  2R  R)  In  R  +  RR  =  -  (p*R  “  Pa 
dR  P  J- 


P  '*  ‘i  r2y_1 


4“  (R2R2  In  R)  =  (2RR2  +  R2  ■— '  )  In  R  +  RR2 
dR  -  qK 


' 


. 
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Substituting  we  get 


2*2 

R  R  In  R 


R 


R, 


R. 


2y 


f  [pAR  -  pi  i  dR  = 


2  p*  ,  2  2. 

p  (VR  ’V 


piRi 


2y 


2~2y 


(R 


2-2y 


\z  2y  )i 


t  = 


R 


R, 


/In~R 


dR 


1 

P 


^1 

Y-l 


r1\2^ 

R~ 


R„ 


^R 


+  p*  51  - 


R, 


R 


T  2 


R 


t  = 


/IrT¥  dR 


R„ 


p*(y-i) 


y2y 

R  , 


V2 

R 


+  1  - 


/  R1 ' 2 
R~ 


(2.7) 


rlR 

For  R  =  R  =  R0  —■  =  0 
max  3  dt 

P* 


Then 


- 

R_1  \2' 

1  - 

R„ 

3  1 

R  \  2y 

R-, 

J, 

1 

K  / 

,  Rn  , 

(y  -  l) 


(2.8) 


Similarly  we  get  for  the  spherical  bubble 

R 


dR 


t  = 


R 


1 


j  Pl 

r 

R-, 

3 

R  \  31  1 

R 

3  1 

\  i 

1 

i  ! 

1 

1 

rH 

! 

>- 

> 

R 

R  J  J 

px 

1  — 

R 

- 

and 


1 


Pi  = 


1 

3 


R 

3 

R. 

1 

1 

u. 

3y 


(y  -  1) 


(2.9) 


(2.10) 


16 


In  either  of  the  two  expressions  for  t,  the  integrals  that 

appear  must  be  evaluated  numerically.  If  the  initial  pressure 

p^  7^  p* ,  the  above  equations  give  the  result  that  the  bubble  will  be 

able  to  oscillate  for  an  infinitely  long  time.  This  is  of  course  not 

true  in  practice  since  the  oscillation  will  be  damped  through  viscous 

losses  in  the  liquid.  For  a  cylindrical  bubble  in  water  with  y  =  1.4, 

=  1  mm  and  =  10  mm  we  get  p^  =  48.6  atm.  and  p^  =  0.075  atm. 

The  graph  of  function  R(t)  for  the  first  half  wave  is  shown  in  Fig. 

(2.4).  The  equation  of  motion  for  a  spherical  bubble  may  be  also 

R  3 

written  in  dimensionless  form..  The  dependance  of  on  t /  (—  pR^)for 
different  values  of  y  and  p^  =  50  atm.  is  shown  in  Figure  (2.5). 

Therefore,  just  the  value  of  the  repetition  frequency  (which 
can  be  easily  measured,  e.g.  piesoelectric  probe)  is  insufficient  for 
the  calculation  of  the  pressure  p.  That  is  why  an  optical  record  is 
necessary. 

A  comparison  between  calculated  values  and  values  obtained 
from  an  optical  record  for  a  spherical  bubble  in  petroleum  after  a 
0.1  joule  discharge  100  ysec  long  is  given  in  Figure  (2.6)  where 
P1  =  40  atm  and  p^  =  0.085  atm.  The  difference  between  the  measured 
and  calculated  values  is  the  result  of  energy  losses  such  as  heat 
exchange,  friction,  vaporization  etc.  Each  energy  loss  increases 
R^  and  decreases  Ry  This  causes  the  change  in  the  time  scale  in 
Figure  (2.5).  Also  any  deviation  from  the  assumed  spherical  shape  will 
result  in  a  difference  between  the  measured  and  calculated  values. 

A  nomogram  has  been  computed  for  the  case  of  a  spherical 
bubble  (Fig.  (2.7))  from  Xvrhich  values  of  p^  may  be  directly  read  for 
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Fig.  (2.4)  The  dependence  of  radius  Fig.  (2.5)  The  dependence  of  R/R  on  t//3/2  pR 

R  on  time  t  for  an  oscillating  for  a  spherical  bubble  in  water/ 

cylindrical  bubble  in  water 
(y  =  1.4) 
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Fig.  (2.7)  The  dependence  of  R„/R^  on 
pressure  p^  for  a  spherical  bubble 
in  water  for  various  values  of  y. 
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different  ratios  of  R./R.  and  different  values  for  y  assuming 

3  4 

p*  =  1  atm.  This  may  be  used  to  obtain  a  quick  overall  picture  by 
evaluating  the  bubble  oscillation  records. 

2.4  The  frictional  losses  of  a  cylindrical  bubble  oscillation 

The  viscous  properties  of  the  liquid  are  also  an  important 
factor  in  the  investigation  of  the  dynamic  processes  of  a  discharge 
between  two  flat  parallel  electrodes.  The  frictional  losses  could  be 
responsible  for  all  of  these  effects: 

(1)  damping  of  oscillations  of  a  simple  bubble  through  the  trans¬ 
formation  of  liquid  kinetic  energy  into  thermal  energy, 

(2)  reduction  of  the  shock  wave  strength,  especially  near  the 
electrode  surface, 

(3)  increased  pressure  in  the  plasma  channel  during  its  expansion 
because  of  the  increase  in  the  force  required  to  move  the 
surrounding  liquid. 

The  exact  description  of  the  behaviour  of  the  viscid  fluid 
around  the  expanding  or  collapsing  discharge  channel  is  a  rather 
complicated  problem.  In  this  chapter  an  approximate  solution  will  be 
obtained,  which  could  give  some  idea  of  the  significance  of  liquid 
viscosity  in  our  problem. 

The  cylindrical  approximation  for  the  oscillating  bubble 
holds  well  in  the  middle  part  of  the  discharge  column.  Near  the 
electrodes  this  approximation  is  probably  inaccurate  even  in  case 
when  liquid  viscosity  is  ignored.  The  theory  of  viscous  flow  gives 
the  frictional  losses  to  be  greatest  near  the  edge  of  the  gas  cylinder, 
where,  if  the  Reynold's  number  is  sufficiently  large,  the  flow  is 
laminar  with  the  appropriate  boundary  layer  with  a  shear  type  of  flow. 
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In  this  case  the  faster  the  flow,  the  thinner  is  the  boundary  layer 
and  the  larger  are  the  friction  losses  per  unit  area.  The  thickness 
of  the  boundary  layer  in  such  a  case  rises  with  increasing  distance 
from  the  middle  of  the  cylinder.  For  distances  from  the  middle  of 
the  cylinder  which  are  larger  than  the  distance  r*  in  which  the 
boundary  larger  thickness  is  equal  to  the  half  of  the  distance  between 
electrodes,  the  flow  is  no  longer  "laminar  with  a  boundary  layer" 
but  instead  "laminar  viscous  flow  with  small  Reynolds  number" 

(Fig.  (2.8)). 


Fig.  (2.8)  The  viscous  flow  around  the  oscillating  cylindrical  bubble 
The  problem  will  be  solved  separately  for 

(a)  r  >  r*  -  case  of  radial  laminar  viscous  flow  between  two  plates 

with  a  small  Reynolds  number". 

(b)  r  <  r*  -  case  of  the  radial  viscous  flow  with  a  boundary  layer. 

The  following  well  known  formula  describes  the  uniform  Laminar 
viscous  flow  with  a  small  Reynolds  number  between  two  plates : 

u(y)  -  %  (d2  -  y2) 

where  2d  -  distance  between  plates 


y  -  viscosity 


. 
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y  -  coordinate  perpendicular  to  plates 
c  -  constant  proportional  to  the  fluid  velocity 
The  average  velocity 


U  = 


_JL 

2d 


2 

3 


U 

max 


U 

max 


For  the  cylindrically  symmetric  case,  in  Fig.  (2.8), 


U(r)  =  UR  f 
Umax  -  !  UR  7 


c  =  U 

2  max 
d 


=  3U„  *  2- 
R  r  ,2 
d 


The  viscous  shearing  stress  tensor 


(2.11) 


13  u.  ,2cy.  ,  OT,  R  y 

t  =  h(y— )  =  uC-v-2-)  =  cd  =  3U  --v 

o  3y  y— d  2p  y=d  ^  ^ 

The  total  frictional  force  for  r  >  r* 

r  _  r 


F  =  2 
a 


2Trrx  dr  = 
o 


12UT)ttRp 

K 


12U  Rrp 

dr  =  — h -  (r-r*)  (2.12) 


r  * 


so  for  r  F  ->  00 

a 

For  case  (b)  where  r  <  r*  we  will  use  the  boundary  layer 

theory.  The  direct  calculation  through  the  evaluation  of  the  boundary 

layer  thickness  is  not  convenience  because  it  leads  to  an  infinite  value 

for  t  in  the  case  where  the  boundary  layer  thickness  5  is  equal  to 
o 

zero  i.e.  for  r  =  R,  A  more  suitable  method  is  the  Karman-Pohlhausen ' s 


momentum  integral  technique  which  allows  the  elimination  of  the 
boundary  layer  thickness  before  integration.  This  method  is  adapted 
below  for  the  case  of  radial  flow. 


. 


- 
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From  Euler’s  momentum  theorem 


5X-  1 


ii  2  2  2  2 

puq  cos  y  ds  F  |  |  u ;  q  =  u  +  v  +  w 


where  £  F  -  the  sum  of  the  external  forces  acting  on  the  fluid  in 
the  considered  volume  q  -  the  velocity  vector,  u,  v,  w,  -components 
of  vector  q. 


Fig.  (2.9)  The  segment  of  the  moving  anulus  and  its  surface  element 

The  forces  acting  on  the  segment  of  the  anulus  in  Figure  (2.9) 
will  now  be  examined. 

The  flow  through  walls 


abef 

cdgh 

bdhg 

aceg 


0  (rigid  wall) 

0  (because  of  symmetry) 


0 

0 


\  g  cos  y  =  0 


(  radial 

flow)  / 

The  total  radial  force  acting  on  the  cylinder 
2tt  2tt 


Fb  =  2  ( 


Fabcd  d*  +  Fefgh  d« 
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The  factor  of  2  accounts  for  the  bottom  and  top  parts  of  the  cylinder. 

d 

f  2 

Ij-  =  2ttR  /  pug  cos  y  dy  =  pU^  2iTRd 
o 

since  u(y)  =  =  constant  at  the  interface  between  the  expanding 

bubble  and  the  liquid. 

Now  we  assume  the  velocity  distribution  in  the  boundary 
layer  to  be 

u  _  .  /IT  y  N 

F - Sln  (2 

max 

then  similarly  as  in  the  previous  case 

6 

u  =  ~  [2U  (d— 6 )  +  2  /  U  sin  £  dy] 

2d  max'  J  max  2  6  J 

o 

Where  the  first  term  in  the  brackets  contains  the  velocity 
profile  in  the  area  between  the  boundary  layers  and  the  second 
contains  the  velocity  profile  within  the  boundary  layer 


also 

Thus 


u 


u 


~  [2U  {d  +  6 (—  -  1)}]  =  U  [1  +  -  1)] 

2d  max  it  '  max'  dm 

v 

-  from  the  continuity  equation. 


U 

max 


UrR 


r[l+|(f-  1)] 


(2.13) 


and  for  r  =  r* 

<5  = 

d 

so  that 

U 

max 

r=r* 

urr 

TT 

urr 

v 

1.57  — r 
r* 

r*i 

TT 

2 

r* 

Equation  (2.11)  gives 
u  R 

U  =1.5  — v 
max  r* 

r=r* 
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The  difference  m  the  values  for  U  is  caused  by  the  assumption 

max  ,  y 

r=r* 

of  a  different  kind  of  velocity  distribution  approximation  in  the 
boundary  layer. 


Then 


I2  =  2irr* 


a 

•/ 


o 


2  2  TT 

pU  sin  (7-  -7)  dy  =  7ir*pU 

max  2d 

r-r* 


max 
r=r* 


d  = 


3  2  2 

ILR  2  ir  i  R  pd 

/TT  R  N  J.  ,  R 

(it*5  1Tr"pd  =  — 75 — 


Fb  -  4,RUR2pd  [||  -  1] 

The  radius  r*  is  dependent  on  the  viscosity  and  the  velocity 

IL,  (the  smaller  the  value  of  r*  the  larger  is  p  and  the  smaller  is  U  ) . 

The  value  of  r*  will  now  be  calculated.  Because  r  ->  r*  if 

5  =  d, the  shape  of  the  boundary  layer  is  required.  As  in  the  case 
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of  a  uniform  flow  past  a  flat  plate  ,  the  following  approximation 
is  made  for  6 . 


6  =  5 


P (r-R) 

P  u 


max 


Then  from  equation  (2.13) 


U  (r) 
max 


urr 


r[l  + 


5(“  "  1) 

TT 


P  (r-R) 


P  U 


max 


or 


U  +  /U 
max  max 


5(4  -  1) 


a  =  b 


a  = 


-  (r-R)  ;  b  = 

P  r 


where 


The  square  of  the  above  equation  gives 

U2  -  (2b  +  a2)  U  +  b2  =  0 
max  max 


(2.14) 
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Equation  (2.14)  also  corresponds  to  the  equation 


U  -  /u  a  =  b 
max  max 

but  for  b  >  0;  a  <  0  and  U  >  0 
—  —  max  — 


the  required'  solution  is 

2 


U  =  b  +  ~~  +  \/a2b  +  ~ 

max  2  v  4 


By  substituting  for  a  and  b 


U 


max 


UR  25p(r-R)(|~l)2 

r  2  ,2 

d  p 


25p  (r-R)  (f-l)2 


j2 

d  p 


U  R  1  25p  (r-R)  C~— 1)  2 

K  ,  -L  IT 

~  +  4  T 

d  P 


But  case  (a)  (laminar  flow)  gives 


U 


UR 
R  7T 


at  r  =  r5 


max  r*  2 

By  substituting  this  condition  into  the  previous  equation  we  finally 
get 


UR 
R  ,  TT 

r*(2 


-  1) 


,  25p  (r*-R)  (~— 1)  2 

i.  _____  TT 

2  ,2 
d  p 


From  the  above  equation  a  value  for  r*  may  be  obtained 

through  the  use  of  a  numerical  iteration  method.  The  force  F  =  +  F 

represents  the  total  radial  frictional  force  which  acts  on  the  cylinder 

of  radius  R,  and  length  2d  expanding  or  contracting  radially  with 

velocity  U  .  The  integral  denoted  by  F^  diverges  as  r  00 .  However, 

because  of  the  slight  compressibility  in  the  liquid,  deformations  and 

finite  dimensions  of  the  electrodes  the  force  F  is  in  fact  finite. 

a 

To  get  a  reasonable  estimate  of  the  friction  force,  the  integration 

will  be  performed  for  r  <  r  only,  where  r  =  1  HO  cm.  For 

a  typical  case  R  =  0.1  cm,  2d  =  0.1  cm,  U^  =  5,000  cm/sec, 

p  =  0.003  cm2/sec.  and  T  =  100°C.  Thus  the  Reynold's  number  is 
H20 


. 
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R  = 


URdp 

y 


5,000  X  0.05 

0.003 


=  83,300 


The  viscous  energy  losses  during  a  half  period  of  oscillation 
will  be  calculated  numerically  under  the  following  assumptions. 

(1)  the  total  amount  of  viscous  losses  during  a  half  period  is  much 
smaller  than  is  the  amount  of  mechanical  energy  stored  in  oscil¬ 
lation  bubble.  Under  this  assumption  the  shape  of  the  oscillation 
bubble  does  not  change  appreciably  and  therefore  calculations  will 

be  made  with  values  U  ,  t,  T  holding  for  inviscid  oscillations. 

K 

(2)  friction  losses  will  be  calculated  for  small  time  intervals  assum¬ 
ing  that  during  this  interval  U  ,  t,  T  are  constant. 

The  viscous  losses  of  an  oscillating  cylindrical  bubble  whose 
radial  time  dependence  is  given  in  Figure  (2.5)  have  been  calculated 
numerically  for  water  and  silicon  oil.  The  application  of  the  Newton- 


Raphson  method  for  the  solution  of  r*  in  equation  (2.15)  gives 


UR 

1.65A(r.*-R)  +  0.57 

i  r . 


+  1.81B 


r »  =  r  ' 

i+1  i 


URR  1  01  9 

1.65A  +0.57  -~2  +  [1.65A +r*-R)  + 


ar2~u 


R- 


r . 
i 


2B 


r*2 


where 


U 


,2 

d  p 


B 


UR 

A(r*-R)[-~ :+  0.825A  (r*-R)  ] 


The  dependence  of  R^r*  and  w  on  t  is  shown  in  Figure  (2.10) 
for  the  case  where  water  and  cylindrical  electrodes  with  a  diameter  of 
10  cm  are  used,  w  represents  the  total  frictional  loss  up  to  time  t, 


that  is 


w 


R 


R„ 


(F  +  Ft  )  dR 
a  b 


where  R  =  R(t) 


. 
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Fig.  (2.10)  The  dependence  of  radii  R  and  r*,  and  frictional 
losses  w  for  an  oscillating  cylindrical  bubble  in  water. 
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The  total  energy  of  the  oscillating  bubble  is  equal  to  half  of  that 
which  is  done  by  an  adiabatic  expansion  of  the  bubble  from  a  diameter 
to  Ry 

WT  =  2  (P3R3  "  P1R1  )  2lTd 

3 

For  the  case  considered  in  Figure  (2.10)  with  2d  =  1mm,  w^  =  1.05  x  10 

dyn  cm  the  frictional  losses  during  a  half  period  of  oscillation  for 

2 

water  and  silicon  oil  (p  =  1.7  g/sec  cm)  are  8.28  x  10  ergs  and 
2.93  x  10^  ergs  respectively.  Thus  the  initial  assumptions  are  invalid 
because  the  energy  losses  are  too  large.  Therefore  the  frictional  forces 
must  be  considered  through  the  numerical  integration  of  the  equation 


(2.7). 

2.5  Dissociation  of  the  water  vapors 

For  many  important  investigations,  for  instance  of  diffusion, 
thermal  conductivity  etc. ,  the  composition  of  the  gasses  contained  in 
the  bubble  is  required.  In  this  section  the  simple  case  of  dissocia¬ 
tion  of  water  vapours  will  be  investigated.  For  simplicity  only  three 
simple  dissociation  reactions  (no  ionization)  will  be  considered 

H20  t  H2+io2 
H2  t  2H 

02  t  20 

Other  reactions  such  as 

OH  t  !°2+K 

O3  —  0  2 


will  be  neglected. 


The  calculation  of  the  degree  of  reaction  e  is  well  described 


in  the  literature29.  The  equilibrium  state  for  the  reaction 
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aA  +  bB  +  cC  +  dD 

is  described  by  the  law  of  mass  action: 


c  d 

XC  XD  c+d-a-b 

Tv  ' 

A  B 


=  K 


where  the  X's  are  mole  fractions  (functions  of  e) 
a,b,c,d  -  stoichiometric  coefficients 
P  -  pressure 

K  -  the  equilibrium  constant 

Values  of  K  are  tabulated  in  some  publications  for  various 
temperatures  and  pressures.  Unfortunately  this  information  is  not 
complete  and  it  is  necessary  to  resort  to  extrapolation.  However, 
some  of  this  data  is  tabulated  for  the  inverse  reactions.  In  such 
cases , values  of  K  for  the  normal  type  of  reaction  may  be  calculated. 
That  is,  if  K  relates  to  the  reaction 
aA  +  bB  cG  +  dD 
and  if  K  corresponds  to  the  reaction 
cC  +  dD  aA  +  bB 

the  formal  relationship  betx^een  K  and  K  is  as  follows : 


K  =  f (e) 


f(l-e) 


If  f(e)  and  K  are  known ;  the  value  of  K  may  be  easily  calculated. 
For  the  reaction 

h2°  t  h2  +  \  02 

the  form  of  the  law  of  mass  action  is  as  follows 


e 

(2+e) 


3/2 

1/2 


(1-e) 


. 


. 

' 
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or 


e (3-e) 


1/2 


(1-E)3/2  p1/2 


=  K 


The  Newton-Raphson  method  will  be  used  for  calculation  of 
values  for  e  from  the  above  equations  on  a  computer 


F(X.) 

Xi+1  =  Xi  “  F' (X.) 

1 


If  a  value  for  K  is  available  for  water,  then  the  use  of 


the  iteration  formula  gives 


:.+1  =  o.  +  /P  [K  -  \j  3  y  •  ]\/(l-e)5(3-e) 


(1-e)  P 


If  the  value  of  K  is  available,  the  iteration  formula  gives 


3/2 


£..,=£.+  [K 
l+l  l 


(2+e) 1/2 (1-e) 


pl/ 2  j  (2+e) 3/2 (1-e) 2 


e1/2[3+  (l-e)]P1/2 


Because  the  graph  of  logpQ  K  against  1/T  (Figure  (2.11))  is 
approximately  a  straight  line,  K  may  be  expressed  as  function  of  T 
as  follows 

log10  K=alTtbl 

where  a^  and  bp  are  constants  and  T  is  in  °K  or 


K  =  exp  [2.30261  x  (—  +  b  )] 

When  P  =  1  atmosphere,  Figure  (2.11)  gives 

3 


a.  =  -13.1  x  10' 


bp  =  +3.03 


Similarly  for  the  reaction 

i  °2  ?  20 2 


a2  =  12.8  x  10' 


b„  =  2.1 


K0  =  / _ 

/l-e2 


2e  pl/ 2 


IC 


.1/2 


/(1-e2)3 


■ 
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,  _  P  ,  IfK  _  2e  1/2  J  (1-e2)3 

£i+l  -  £1  +  2  K2  ^  )  V  P 


and  for  the  reaction 


H2  t  2H 


K3 


4e2P 

1-e2 


a3  =  -23.48  x  10' 


=  6.35 


K3  = 


8e 


n  2.2 

(1-e  ) 


2  2  2 
e  =  e  +  (k  -  AcJP)  OzO_ 

1+1  1  U3  ,  2’  8e 

1-e 


From  the  above  relations,  values  for  e 
be  calculated  for  particular  values  of  the  temperature  T  and  the 


H20’  £H2’  £02/2  may 


pressure  P.  Finally,  the  complete  reaction  is  described  by  the 
following  equation, 


n, H  O  noHo0  +  n0H0  +  n,H  +  nc;0o  +  n.O 

12  22  32  4  52  6 


(2.15) 


The  reaction  described  by  equation  (2.15)  is  depicted  graphically 
below 


nl('1  £H2(PH2° 


n  e  H  +  —  n  e  0 
1  H20  2  2  FH  O  2 


nleH20('1  £H2^H2  2nl£H20£H2H  2  nl£H20(1  £02^°2  nl£H20£02° 


The  dependance  of  the  equilibrium  constant  e  on  the  temperature  for  the 

three  basic  reactions,  with  P  =  1  atm.,  is  shown  in  the  Figure  (2.12). 

2.6  Ionization  and  thermal  conductivity  of  the  water  vapours 

The  method  for  calculating  the  degree  of  ionization  as  well 

as  the  thermal  conductivity  of  the  water  vapours  is  described  in 

20 

detail  by  Gross  et  ah'  .  For  a  complicated  mixture  of  different 

20 

kinds  of  molecules,  atoms,  ions  and  electrons.  Gross  et  al .  derive 
from  the  system  of  Saha’s  equations  and  equations  of  state  (under 
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Fig.  (2.11)  The  dependence  of  log^K  on  1/T  for  reactions 

(A)  H20  t  H2  +  02/2 

(B)  H2  t  2H 

(C)  0 2 / 2  t  0 
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Fig.  (2.12)  The  dependence  of  equilibrium  constant  e  on  temperature  T  for 
the  reactions  (A)  l^O  1-  H2  +  O2/2  (B)  H2  £  aH  (C)  O2/2  Z  0  at 

1  atmosphere  pressure. 
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some  simplif icating  assumptions)  the  following  expressions  for  the 
particle  concentrations: 


VT) 

ne  kl' 


[A  + 


VT) 


n 


e  2kT 


V(T) 

p  3+8  — - 

3+7  — " — 


-  1]  for  T  <  2  x  10  °K 


for  T  >  2  x  104oK 


n 


S0i-+(T) 
P  4+9  -~— 


e  kT 


1+2 


s0-h-(t) 


for  T  >  4.5  x  10  °K 


2  P  2SH(T) 


nH  3  kT  ^ 1  P 


/l  ■  P 


(/I  + 


sh(t) 


"  D] 


2  1  "  r/l  +  -  1] 


nH+  3  kT  SH(T)  [  1  +  Sr(T) 


P  So(T) 

n  =  [1  -  2 


o  3kT 


P  ^  +  S  (T) 
o 


3+8 


V(I> 


n 


o+  3kT 


S0+(T)  S  +(I)  2 
18+90  — -  +98 (-2= — ) 


So+(T)  So+(T)  2 

p  6  -2-g-  +  14(-V+> 


rlo++  3kT 


So+(T)  V(T)  2 

18+90  -2= - +  98  ( — 5 - ) 


+  5(- 


o 


++ 


(T)  2 


no++-f  3kT 


So++(T) 

4+16  -+5 - f-  14  ( 


S0++(T)  2 

A - ) 


where  P  -  pressure;  ng,  n^,  nQ,  nQ+>  n0-H->  %+++  “  concentrations 

of  electrons,  hydrogen  and  oxygen  atoms  and  ions  respectively;  S^,  S^-f, 
and  Sq++  -  right  hand  sides  of  the  Saha  equations,  which  have  the 
general  form 
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ni+l 


2Z .  (27Tm  ) 

l+l  e 


3/2 


e . 
l 


n.  Pe  ”  Z. 


(kT)5/2  e  kT 


where  n.  ,  n .  ,  ,  -  concentrations  of  particles  in  the  i'th  and  i+l'th 
l  l+l 


degree  of  ionization;  Z Z._^  -  partition  functions  of  particles  in 

i'th  and  i+l'th  degree  of  ionization,  pg  -  partial  pressure  of 

electrons,  m  -  mass  of  the  electron,  h  -  Planck's  constants,  k  - 
e 

Boltzman's  constant,  T  -  temperature,  e.  -  ionization  energy  for  the 
i'th  degree  of  ionization. 

The  expressions  for  n  and  n^(T  >4.5  x  10  °K)  have 

been  derived  by  the  author  by  using  a  method  similar  to  that  used 

,  „  ,20 

by  Gross  et  al 

The  thermal  conductivity  of  the  gas  inside  the  bubble  is 
now  determined  for  a  ten  component  mixture  with  the  concentrations 

nH20’  nH2’  n02’  V  V  nH+’  n0+'  V^’  V++’  "e  aS  funcClons  o£  the 
pressure  p  and  temperature  T.  The  method  for  calculating  the  thermal 

20 

conductivity  is  again  described  by  Gross  et  al  where  the  expression 

for  the  thermal  conductivity  is  given  by 

k 


KN  = 


n .  r 

~~  c.  A.  (3r  C.  +  C.  .  J 

3  i  i  k2  itrans  i  mt 


i=l 


The  thermal  conductivity  associated  with  the  diffusion  of  dissociation 
(the  diffusion  of  ionization  is  neglected)  kd  is  given  by 


K  = 
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h  .  p  .  „  , 

JU  [Mi  (m 

A  L  ST 
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U  1 

m.  — )  -  — 
A  e2  n. 


3p- 
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h  .  p  .  Sp 
1  1  [ 
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where  k  -  number  of  assumed  components;  c^  -  mean  thermal  velocity, 

H.  -  mean  free  path,  C.  .  and  C.  .  -  specific  heats  of  trans- 

i  v  *  1  trans  1  mt 

itional  and  internal  degrees  of  freedom  for  the  i’th  component, 

h  -  enthalpy  of  the  i’th  component,  p.  -  partial  density  of  the  i ' th 

i  1 

component,  A  =  n^e^  +  n^^*  ^  =  nMe2  nAe3’  anc^exes  ^  an<^  ^  re^er 

to  molecules  and  atoms  respectively.  The  factors  e^,  e ^ represent 
the  mutual  collisions  of  molecules,  atoms  and  atom-molecule  collisions 

respectively . 

The  total  thermal  conductivity  k  is  given  by 

k  =  k  +  k 

N  D 

Simplifying  assumptions  such  as  the  neglect  of  the  higher 
orders  of  ionization  at  low  temperatures  and  the  neglect  of  the  in¬ 
fluence  of  molecules  and  atoms  at  higher  temperatures  were  made  in 
the  computation  of  k. 

2.7  The  temperature  distribution  and  heat  outflow  in  the  space  between 
the  plasma  and  the  liquid 

By  assuming  that  the  plasma  column  is  described  by  the  theory 

1  7 

of  a  thermally  inhomogeneous  plasma  column  (see  also  section  2.2), 
the  temperature  distribution  in  the  space  between  plasma  and  liquid 
can  be  calculated.  The  temperature  gradient  in  a  cylindrical  shell 
of  a  length  1  is: 

dT  Q 

dr  2tt1  k 

c 

where  Q  -  constant  rate  of  heat  outflow,  T  -  temperature,  r  -  radius, 
k  -  thermal  conductivity. 

The  direct  integration  of  the  above  equation  in  the  space 
between  liquid  and  plasma  is  not  possible  because  k  is  a  complicated 
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function  of  T  (assuming  the  pressure  P  to  be  a  constant) .  A  numerical 
iteration  technique  must  therefore  be  used  to  solve  the  boundary  value 
problem  defined  by  the  temperature  of  the  walls  of  the  bubble  and  the 

plasma  column. 

The  temperature  distribution,  the  thermal  conductivity  and 
Q  have  been  computed  in  Figure  (2.13)  for  an  approximation  to  the 

state  existing  at  the  end  of  a  100  ysec  discharge  shown  in  Figure 

•  8 

(3.1a).  For  this  case  Q  =  1.9  x  10  erg/sec.  This  can  be  compared 

with  the  energy  of  the  oscillating  bubble  described  in  section  2.5 

3 

which  has  a  kinetic  energy  of  1.05  x  10  erg  and  an  oscillation  period 
of  3.14  msec.  The  energy  outflow  during  the  first  microsecond 
(assuming  that  temperature  inside  the  bubble  does  not  change)  is  com¬ 
parable  with  the  kinetic  energy  of  the  bubble.  However,  the  temp¬ 
erature  and  therefore  the  heat  outflow  drop  rapidly  after  the  dis¬ 
charge  extinguishes.  Thus  the  exact  solution  of  the  behaviour  of 
the  cylindrical  or  spherical  bubble  after  the  end  of  the  discharge 
requires  the  consideration  of  all  thermal  effects.  This  could  be 
obtained  by  a  numerical  solution  for  the  nonlinear  partial  differ¬ 


ential  equation. 
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CHAPTER  3 

INSTABILITY  OF  THE  PLASMA  COLUMN 

3.1  Unstable  electrode  phenomena 

The  presence  of  an  electrode  surface  causes  a  cathodic 

sheath  to  appear  at  one  end  of  the  discharge  column  and  a  anodic 

sheath  to  appear  at  the  other  end.  Each  sheath  has  its  mechanism 

characterized  by  a  higher  electric  field  intensity  than  that  found 

in  the  remainder  of  the  column  and  could  under  certain  conditions 

influence  the  behaviour  of  the  discharge  column. 

The  motion  and  division  of  the  cathodic  spot  has  been 

observed  for  some  time.  However,  the  theoretical  description  of  this 

phenomena  has  not  been  sufficiently  developed.  Because  in  many  spark 

12  2 

discharges  the  current  density  reaches  values  of  up  to  10  amps /mm  , 

some  cathodic  spot  emission  mechanism  is  probable.  The  extremely 
high  emission  current  density  in  the  cathodic  spot  cannot  in  all 
cases  be  explained  by  thermic  or  field  emission. 

The  Longini's  model^  best  explains  the  mechanism  of  cathodic 
spot  emission.  According  to  the  Longini’s  model  the  very  fast  electrons 
disturb  the  upper  layers  of  the  cathode  ion  lattice  in  a  fashion  similar 
to  that  by  ambipolar  diffusion.  A  special  transition  layer  is  created 
by  the  cathodic  material  leaving  the  cathode  with  a  metallic  kind  of 
electrical  conductivity  and  a  very  low  emission  work  function.  Accord¬ 
ing  to  Steenbeck's  minimum  principle  the  discharge  current  is  concentrated 
at  the  cathodic  spots.  By  using  this  principle  it  is  also  possible  to 
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deduce  that  the  number  of  cathodic  spots  is  proportional  to  the  current. 
The  cathodic  spots  have  been  observed  to  move  more  or  less  at  random 
on  the  surface  of  the  cathode.  The  explanation  of  this  motion  is  un¬ 
fortunately  unknown.  A  remarkable  feature  is  that  the  stagnation  of 
the  spot  leads  to  a  slight  increase  in  the  cathodic  voltage  drop. 

In  some  types  of  arcs,  the  current  density  is  very  unevenly 
distributed  on  the  surface  of  the  anode.  This  results  in  the  current 
being  collected  by  anodic  spots  which  may  sometimes  be  indistinguishable 
from  cathodic  spots. 

The  seige  of  the  relaxation  phenomena  in  the  anodic  region 
is  commonly  termed  as  "hash"  by  engineers.  It  occurs  if  the  voltage 
drop  in  the  anodic  region  is  large.  The  anodic  sheath  is  virtually 
short  circuited  by  a  plasma  bridge,  which  is  created  by  bursts  of 
ionization . 

3.2  Unstable  column  phenomena 

Theoretical  studies  of  the  macroscopic  instability  of  a 

plasma  column  in  the  presence  of  an  external  magnetic  field  have  been 

21  22  23 

made  by  Kadomtsev  and  Nedospasov  '  ,  Hoh~~,  Paulikes,  Pyle  ,  and 

r\  / 

Schmidt"  .  Instabilities  similar  to  those  studied  in  the  case  of  an 

2 

external  magnetic  field  have  been  observed  by  Kocian  in  the  field 
free  case.  Similar  instabilities  were  observed  through  the  use  of  a 
time  resolving  camera  by  the  author  in  discharges  in  a  liquid,  Figure 

3.1b. 

The  instabilities  of  self  pinched  discharges  are  described 

q  r  26  23  28 

by  Cruzon"  ,  Wyld  and  Kadomtsev  .  Kadomtsev  gives  the  following 

criterion  for  the  onset  of  the  pinch  instability  I 
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Fig.  (3.1)  Records  of  discharges  in  liquids  made  by  a  time  resolving  camera 
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d  In  r  2+88 

2 

where  3  =  8rrp/B  is  the  ratio  of  the  plasma  pressure  to  the  magnetic 
pressure  and  the  other  symbols  as  described  in  the  previous  chapters. 

The  above  conditions  imply  that  for  the  pinch  instability 
not  to  occur,  the  plasma  pressure  must  not  diminish  too  rapidly  with 
increasing  r  (a  condition  well  satisfied  with  the  discharges  under 
consideration).  If  the  pressure  falls  off  more  rapidly,  the  pinch 
is  unstable  to  an  axially  symmetric  perturbation  and  a  "necking"  or 
"sausage"  instability  will  occur. 

28 

Thermal  instabilities  are  described  by  Kadomtsev  and  by 

18  ?  8 

Ecker  .  Kadomtsev  describes  the  case  with  an  external  magnetic 
field,  while  Ecker  describes  the  field  free  case.  In  addition,  Ecker 
describes  the  influence  of  the  outer  electrical  circuitry  on  the  in¬ 
stabilities.  The  essence  of  his  theory  which  is  based  on  the  thermally 
inhomogeneous  plasma  column  model^  (see  equations  (2.4),  (2.5),  (2.6)) 
is  as  follows.  A  density  disturbance  produces  a  temperature  disturbance. 
The  temperature  disturbance  in  turn  results  in  a  change  of  the  ion 
production  coefficient  so  that  a  coupling  between  the  density  and 
temperature  exists.  Secondly  a  density  disturbance  will  cause  a  change 
in  the  longitudinal  electric  field.  This  too  will  influence  the  in¬ 
stability  via  the  particle  production  and  particle  losses.  The  influence 
of  the  outer  circuit  is  taken  into  account  by  the  simple  relation 

U  =  I  (ft  +  n.)  =  Ifi  +  Elc  (3.2) 

o  1  o 

where  U  -  constant  voltage  supply,  -  outer  resistance,  fi.  -  resist¬ 
ance  of  the  discharge,  1c  -  column  length  and  E  -  electrical  field. 
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By  using  a  first  order  theory  to  describe  the  perturbations 

in  density  and  temperature  which  are  assumed  to  propagate  at  a  subsonic 

18 

velocity,  Ecker  obtains  the  following  criterion  for  the  onset  of  a 
thermal  instability. 
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where  I __  -  critical  current,  X  and  T  are  heat  conductivity  and 
cr  oo  oo  J 

temperature  within  the  Shottky approximation  respectively. 

Without  an  external  resistance  (a-*o)  the  plasma  column  is 
always  unstable.  With  increasing  resistance  (a->l)  the  column  is 
stabilized  for  currents  below  the  critical  value  I 


cr 


3.3  Other  unstable  phenomena 

The  exploding  liquid  and  the  vapourized  metals  streaming  from 
the  electrode  surface  can  build  up  additional  branches  of  the  main 
channel.  Traces  of  such  branches  have  been  observed  on  the  polished 
electrode  surface  after  a  discharge  in  a  liquid. 

The  vaporization  of  the  electrode  surface  coupled  with  the 
irregularity  of  the  moving  liquid  can  cause  the  sudden  changes  of 
the  pressure  inside  the  gas  bubble.  A  rough  and  irregular  electrode 
surface  in  the  proximity  of  the  discharge  can  strengthen  the  above 


effect. 
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3.4  Correction  of  the  thermal  instability  condition 

From  the  thermal  instability  theory  the  critical  current 
Icr  is  defined  by  equation  (3.4).  In  this  expression  a  is  assumed 
to  be  a  constant.  According  to  the  formula  (3.5)  a  is  a  constant  as 
far  as  the  resistance  of  the  discharge  is  constant.  But  from 
previous  experiments,  is  known  that  the  voltage  on  the  spark  gap  is 
approximately  constant  for  a  very  wide  range  of  discharge  currents. 
This  will  be  further  discussed  in  chapter  4.  Thus  the  channel 
resistivity  CL  is  dependent  on  the  current  I.  The  following  approxi¬ 
mate  relation  for  C  therefore  holds: 
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where  3  is  a  constant  and  from  equation  (3.5) 
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(3.9) 


The  cases  where  a  =  constant  and  where  a  ^  constant  are 
illustrated  in  Figure  3.2a  and  3.2b  respectively.  In  the  case  where 
a  =  constant,  a  decrease  in  the  current  I  results  in  a  transition 
from  the  unstable  region  to  the  stable  region  (Fig.  3.2a).  If  the 
stabilizing  character  of  the  outer  circuitry  is  preserved  in  the  case 
where  a  t  constant,  a  decrease  in  the  current  can  lead  to  a  transi¬ 
tion  from  the  stable  to  the  unstable  region  (Fig.  3.2b). 
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Fig.  (3.2)  The  dependence  of  I 

(a)  a  =  const 

(b)  a  4-  const 


on  I  for 
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As  in  the  work  by  Ecker  '  ,  the  coefficient  A  in  equations 
(3.4)  and  (3.9)  is  assumed  constant,  an  assumption  that  does  not 
always  hold  true.  However,  the  object  of  the  analysis  is  to  show 
that  indeed  a  transition  from  the  stable  to  an  unstable  region  is 
possible  by  decreasing  the  current  I.  The  above  conditions  may  be 
similarly  applied  to  the  higher  modes  of  the  instability. 

3.5  Summary 

From  the  previous  chapters  it  follows  that  there  are  many 
possible  reasons  for  the  occurrence  of  discharge  instabilities  in 
liquids.  The  theoretical  analysis  of  these  instabilities  is  a 
difficult  problem.  In  many  cases  not  only  are  the  instabilities 
not  theoretically  described  as  yet  but  the  exact  physical  explana¬ 
tion  of  the  discussed  phenomena  is  not  known,  for  instance  the 
creation  of  a  cathodic  spot  and  its  motion. 

A  more  detailed  experimental  analysis  is  possible  however. 
Figure  3.1  contains  some  records  made  by  a  time  resolving  camera  of 
discharges  in  a  liquid.  The  bunching  which  appears  in  Figure  3.1b 
is  similar  to  a  pinched  column.  The  origin  of  the  other  irregularities 
may  be  due  to  thermal  instabilities.  The  middle  part  of  Figure  3.1b 
is  similar  to  the  "schlauch"  instability  (or  a  helix  instability) 
described  by  Kocian  for  a  low  pressure  plasma  column.  The  influence 
of  the  electrode  phenomena  is  not  clear .  Unfortunately  these  records 
are  not  reproducible  so  that  a  more  detailed  analysis  is  impossible. 

Further  conclusions  regarding  discharge  instabilities  will 
be  made  after  the  description  of  the  experiments  which  were  performed 
to  check  the  influence  of  the  outer  electrical  circuitry,  the  geometric 
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arrangraent  and  the  liquid  and  electrode  materials.  For  instance,  in 
conjunction  with  the  thermal  instability  theory,  the  effects  of  both 
the  outer  resistance  damping  and  electrode  phenomena  will  be  invest¬ 


igated. 
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CHAPTER  4 

EXPERIMENTAL  RESULTS 

4.1  Apparatus  description  and  method  of  measurement 

The  spark  gap  consisted  of  two  electrodes  held  by  movable 
copper  electrode  holders,  all  of  which  were  placed  into  a  colorimetric 
cuvette  filled  with  a  liquid.  A  copper  vessel  with  observation  ports 
was  placed  into  the  cuvette  to  protect  the  cuvette  against  the  strong 
shock  waves  created  by  the  discharge  in  the  spark  gap.  The  spark  gap 
area  was  illuminated  from  the  rear  side  by  a  simple  light  source  and 
was  observed  from  the  front  side  through  a  stereoscopic  microscope. 

The  microscopic  ocular  was  equipped  with  a  scale  which  allowed  the 
direct  measurement  of  the  distance  between  the  electrodes.  The  small¬ 
est  measurable  distance  was  25  pm. 

The  spark  gap  was  driven  by  a  simple  pulse  independent  EDM 
generator  which  is  described  in  detail  in  Appendix  A.  An  SCR  connects 
the  spark  gap  with  a  charged  condenser  through  a  load  resistor.  By 
making  the  time  constant  of  this  RC  element  sufficiently  large,  the 
current  streaming  through  the  spark  gap  could  be  assumed  to  be  constant 
during  a  small  portion  of  the  discharge  time.  In  the  case  considered 
the  time  constant  x  =  2.7^70  msec  and  the  observation  time  20  psec. 

The  inductance  of  the  connection  between  the  spark  gap  and 
the  discharge  condensor  was  minimized  by  using  short  thick  wires  in 
a  parallel  line  arrangement.  The  total  inductance  of  the  discharge 
circuit  was  0.713  pH  thus  giving  a  characteristic  frequency  of  12.7  kHz 


for  the  RLC  circuit. 
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SCR’s  with  a  maximum  turn  on  time  of  3  psec  were  sufficient 
for  this  experiment  since  the  object  of  this  study  does  not  include 
the  very  fast  phenomenon  occurring  at  the  outset  and  at  the  end  of  the 
discharge. 

The  discharge  current  and  voltage  were  measured  with  a 
Tektronix  502  oscilloscope.  The  voltage  on  the  spark  gap  was  measured 
directly  by  the  scope  probe.  The  discharge  current  was  measured  on  a 
shunt  which  was  switched  into  the  discharge  circuit.  The  shunt  com¬ 
prised  of  a  single  wire  with  a  resistance  of  0.0775 and  an  inductance 

__  O 

of  2.1  x  10  Henrys.  The  shunt  inductance  was  compensated  for  by  the 
capacitor  bridge  in  the  probe  used  for  measuring  the  shunt  voltage. 

All  electrical  measurements  on  the  discharge  were  made  in  the  time 
interval  between  60  and  80  psec  after  the  onset  of  the  discharge.  The 
maximum  voltage  to  which  the  discharge  capacitor  could  be  charged  was 
Ximited  by  the  forward  breakdown  voltage  of  the  SCR  s.  Foi  the  range 
of  electrode  distances  that  could  be  obtained,  this  voltage  was  in¬ 
sufficient  for  a  pure  liquid  breakdown.  This  problem  was  solved  by 
means  of  a  tiny  electrically  conductive  bridge  which  was  built  up 
between  the  two  electrodes  before  the  discharge.  This  bridge  was 
created  from  a  very  small  amount  of  graphite  powder  which  was  dis¬ 
persed  in  the  liquid.  The  single  graphite  particles  were  attracted  by 
a  weak  electrical  field  into  the  maximal  gradient  space,  where  being 
polarized,  attracted  one  other.  The  diameter  of  a  graphite  bridge 
formed  in  this  way  does  not  exceed  few  microns  and  its  further  influ¬ 
ence  is  neglected. 

In  addition  to  the  discharge  current  and  voltage  amplitudes 


and  waveforms,  the  electrode  distance,  the  value  of  the  loading  resistor 
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and  the  source  voltage  were  recorded.  The  evaluation  of  the 
current  and  voltage  measurements  was  simplified  by  noting  only  the  dc 
component  and  the  average  ac  component  amplitude  in  the  interval 
6CK80  ysec  after  the  onset  of  the  discharge. 

Because  of  a  large  deviation  in  the  values,  each  reading  was 
repeated  about  10  times  and  the  average  value  calculated.  All  experi¬ 
ments  were  performed  at  room  temperature  (^20°C)  and  760  mm  Hg  pressure. 
4.2  Liquid  and  electrode  materials  used 

The  number  of  required  measurements  grows  rapidly  with  the 
number  of  liquids  and  electrode  materials  examined.  In  addition,  the 
electrode  distance,  the  discharge  current  and  load  resistor  can  be 
varied.  Therefore  the  present  investigation  will  only  consider  materials 
and  parameters  which  may  be  readily  obtained  and  described. 

The  following  representatives  of  the  main  liquid  groups 
were  used  in  the  experiments  that  were  performed: 
inorganic:  water  and  nitrogen 

silicate:  silicon  oil 

organic:  benzene  (cyclic  hydrocarbon)  and 

hexane  (saturated  hydrocarbon) 

The  most  detailed  measurements  were  made  with  water  because 
water  is  a  well  defined  substance,  hamless,  cheap  and  is  relatively 
simple  chemically.  However,  the  conductivity  of  water  is  strongly 
enhanced  even  with  a  small  amount  of  impurities  such  as  CO^  or  NaCl. 

This  results  in  a  strong  electrolytic  effect  during  the  formation  of 
the  graphite  bridge.  For  this  reason  the  water  which  was  used  was 
distilled, deionized  and  had  the  CO^  removed  by  boiling.  All  metallic 
parts  immersed  in  water  were  first  passivated  in  HNO^. 
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Liquid  nitrogen  was  also  used  because  of  its  very  low  temp¬ 
erature  and  chemical  simplicity.  The  experiment  was  somewhat  more 
complicated  than  that  in  which  water  was  used  because  of  the  low 
temperature  of  the  nitrogen.  This  necessitated  the  use  of  a  metallic 
container  with  double  observation  windows  as  well  as  continuous  com¬ 
pensation  for  the  expansion  of  the  electrode  holders. 

The  chemically  simple  monoatomic  inorganic  liquids  such  as 
bromine  were  not  used  because  of  their  poisonous  nature. 

Of  the  many  types  of  silicon  oils  available,  the  use  of  the 
Dow  Corning  705  diffusion  pump  oil  (chemically  a  pentaphenyltrimethyl- 
trisiloxane)  was  investigated.  In  addition  to  the  presence  of  in¬ 
organic  silicon  in  chemical  bonds,  a  remarkable  property  of  this 
liquid  is  a  viscosity  independent  of  temperature,  over  a  wide  range. 

Benzene  C^H^.  and  hexane  CH^  (CH^) ^CH^  were  used  as  the 
representatives  of  the  cyclic  and  saturated  hydrocarbons  respectively. 
Mixtures  of  various  liquid  hydrocarbons  are  usually  used  in  EDM. 

Electrodes  made  from  both  copper  and  tungsten  were  tried. 
Differences  in  properties  such  as  melting  point,  emission  potential 
and  thermal  conductivity  between  these  two  materials  allows  an  assess¬ 
ment  of  their  influence  on  unstable  electrode  phenomenon. 

The  shape  of  the  electrodes  used  was  either  that  of  a  flat 
plate  or  a  tip.  The  combination  of  two  flat  electrodes  (corners  on 
left  electrode  are  rounded  to  eliminate  edge  effects)  is  shown  in 
Figure  4.1a.  A  tip  electrode  is  shown  in  Figure  4.1b. 
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Fig.  (4.1)  The  shapes  of  electrodes 
4.3  Experimental  data 

As  discussed  in  chapter  3,  the  experimental  investigation 
is  limited  to  a  characteristic  case  of  interest  to  EDM.  The  flat 
Cu-Cu  electrode  combination  in  water  was  chosen  as  a  basic  configura¬ 
tion  for  reasons  outlined  in  chapter  3. 

Most  of  the  experimental  measurements  were  made  on  the  basic 
configuration  in  which  the  electrode  separation  d,  the  source  voltage 
U  and  the  load  resistor  were  varied.  For  each  case,  values  of 
the  channel  voltages  (U^  and  U^)  and  the  discharge  currents  (I^c 

and  I  ) ,  as  they  are  defined  in  a  typical  record  (Fig.  4.2)  were 
ac 
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With  a  load  resistor  larger  than  200  W,  a  wide  fluctuation 

in  the  measured  values  for  I  was  obtained.  For  this  reason,  the 

ac 

present  analysis  is  restricted  to  smaller  values  for  R  .  The  typical 

J-j 

dependence  of  I  on  I,  for  R  =  36  W  is  depicted  in  Figure  (4.3). 

a  c  dc  L 

The  influence  of  the  electrode  distance  (for  various  values  of  R^)  on 
I  is  slight.  For  this  reason,  the  average  measured  value  of  I 

3.C  clC 

(with  the  electrode  separation  in  the  range  25-100  p)  as  a  function 
of  R^  is  given  in  Figure  (4.4).  The  stabilizing  effect  of  the  load 
resistor  R^  as  well  as  decrease  of  1^  with  an  increasing  1^^  is  shown 
in  Figure  (4.4).  The  irregularity  of  the  curve  for  R^  =  155  W  is 
probably  the  result  of  the  large  deviations  in  the  measured  values. 

Measurements  with  values  of  1^  below  3  amperes  were  not 
made  because  of  discharge  interruptions  as  well  as  breakdown  often  not 
occurring  at  all. 

The  dependence  of  the  voltage  on  the  current  I^c  is 

shown  in  Figure  (4.5).  From  this  figure,  it  may  be  seen  that  the 

assumption  made  in  section  (3.4)  that  the  voltage  is  constant  over 

a  wide  range  of  values  for  the  discharge  current  I  is  valid. 

The  voltage  U,  decreases  with  decreasing  electrode  distance, 
dc 

The  dependence  of  the  average  values  o  f  Udc  on  the  electrode  distance 

is  shown  in  Figure  (4.6).  These  curves  are  approximately  straight 

lines  with  only  a  small  dependence  on  R^.  By  using  the  method  of 

minimal  squares,  the  average  dependence  of  U^c  on  d  is  given  in 

Figure  (4.6)  (dotted  line).  The  extrapolation  of  this  line  gives 

U  =19.06  volts  for  d=0.  This  voltage  represents  the  sum  of  the 
dc 
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die  and  cathodic  voltage  drops.  The  slope  of  the  average  (dotted) 
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Fig.  (4.3)  The  dependence  of  I  on  Id  for  different  electrode 
distances.  The  liquid  dielectric  is  water,  =  360 
and  two  flat  Cu  electrodes  are  used. 
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Fig.  (4.5)  The  dependence  of  on  Idc  for  various  values  of  d. 

All  other  parameters  are  as  given  in  Fig.  4.3. 
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Fig.  (4.6)  The  dependence  of  LTdc  on  the  electrode  distance  d  for 
different  R  .  All  other  parameters  are  as  given  in  Fig.  4.3. 
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line  gives  the  electric  field  along  the  middle  part  of  the  discharge 
column,  which  in  this  case  is  E  =  394.4  volts/cm.  Because  the  slope 
of  the  measured  curves  does  not  change  appreciably  in  the  proximity 
of  d  =  25  y  (and  they  in  fact  must  end  in  origin) ,  the  sum  of  cathodic 
and  anodic  sheath  thicknesses  must  be  less  than  25  p.  (Measurements 
previously  made  by  the  author  suggest  a  sheath  thickness  of  about  5  y) . 
Thus  the  major  portion  of  this  discharge  column  is  a  axially  homogeneous 
cylinder. 

The  dependence  of  I  on  I,  (where  R  =  19  and  the  flat  Cu 

cLC  dc 

electrode  distance  d  —  50  y)  for  various  liquids  is  shown  in  Figure 
(4.7).  Of  the  liquids  examined,  the  most  unstable  discharge  occurred 
in  water.  The  amplitude  of  I  is  smallest  in  the  case  of  a  discharge 

in  silicon  oil  and  for  currents  greater  than  30  amperes  was  less 

than  the  resolution  of  the  measurement.  The  voltage  characteristics 
of  the  liquids  tested  were  practically  identical  with  those  of  x^ater 
in  Figure  (4.5)  . 

The  results  of  measurements  with  four  different  combinations 
of  flat  W  and  Cu  electrodes  (R^  =  19  ft,  d  =  25  y)  in  water  are  shown  in 
Figure  (4.8).  In  general,  the  stability  of  the  discharge  is  reduced 
if  one  of  the  electrodes  is  made  of  tungsten.  The  discharge  is  less 
stable  with  the  Cu_W+  combination  than  with  the  Cu+W_  combination. 

Of  the  four  combinations,  the  WW  configuration  is  the  most  unstable. 

The  results  of  measurements  with  various  combinations  of 
shapes  for  the  Cu  electrodes  (R^  =  19  ft,  d  =  50  y)  in  water  aie  given 
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Fig.  (4.8)  The  dependence  of  I  on  1^  for  various  Cu-W  electrode  combinations. 

All  other  parameters  are  as  given  in  Fig.  4.3. 
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cathode 

anode 

I  [ amp ] 

ac 

flat 

flat 

0.193 

flat 

tip 

0.183 

tip 

tip 

0.416 

tip 

flat 

0.447 

The  increase  of  instabilities  with  a  tip  cathode,  regardless 
of  the  shape  of  anode  is  evident.  The  slight  difference  between  first 
and  second  combination  could  be  due  to  experimental  error. 
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CHAPTER  5 

DISCUSSION  AND  SUMMARY 

5.1  Discussion  of  theoretical  investigation 

One  of  the  most  important  phenomena  in  EDM  is  the  pressure 
in  the  discharge  area.  For  this  reason,  a  theoretical  investigation 
has  been  made  of  the  gas  bubble  oscillations  after  the  discharge 
extinguishes.  If  the  equation  of  motion  together  with  the  damping 
forces  is  known,  the  pressure  which  prevailed  at  the  end  of  the  dis¬ 
charge  inside  the  bubble  may  be  calculated  from  an  optical  record. 
Starting  with  a  long  impulse  which  is  successively  shortened,  the 
pressure  at  the  end  of  each  discharge  can  be  calculated.  In  this 
fashion  the  pressure  profile  for  most  of  the  length  of  discharge  can 
be  obtained. 

In  this  thesis,  the  equation  of  motion  of  an  ideal  gas  bubble 
oscillation  has  been  derived.  The  agreement  of  such  a  description 
with  experimental  data  is  very  good  during  the  second  half  wave. 

(See  Fig.  2.6  and  3.1a)  The  method  for  calculating  frictional  thermal 
energy  losses  has  been  extended.  Thus  a  more  accurate  description  of 
the  bubble  behaviour  immediately  after  the  discharge  extinguishes  is 
presented.  A  more  detailed  investigation  of  the  actual  shape  of  the 
gas  bubble,  which  in  fact  is  not  an  exact  cylinder  or  sphere  would 
further  improve  the  accuracy  of  this  method  for  calculating  the  pressure 
profile. 

However,  even  if  damping  effects  are  neglected,  the  approximate 
calculation  of  the  pressures  inside  the  spark  is  still  very  valuable 
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because  by  noting  that  the  pressure  at  the  end  of  the  discharge  is 
relatively  low  (in  Fig.  (3.1a)  it  is  0.32  atm)  and  that  the  subsequent 
drop  in  pressure  is  gradual  (minimum  in  Fig.  (3.1a)  is  0.085  atm),  it 
may  be  concluded  that  the  transport  of  material  melted  by  the  boiling 
mechanism  occurs  mainly  during  the  discharge  and  not  after  the  discharge 
extinguishes.  This  is  in  contrast  to  the  mechanism  occurring  with 
relaxation  generator  EDM  process. 

5.2  Discussion  of  the  experimental  data  and  comparison  with  theory 

The  experimental  data  in  chapter  4  represents  over  4,000 
measurements.  The  poor  reproducibility  of  individual  measurements 
necessitated  a  large  number  of  measurements  followed  by  a  statistical 
evaluation.  Even  if  the  absolute  values  in  this  data  are  not  consider¬ 
ed  as  exact,  the  general  behaviour  of  discharges  in  liquids  are  never¬ 
theless  still  characterized. 

The  stabilizing  effect  of  the  load  resistor  is  evident  in 
Figure  (4.4)  and  in  good  agreement  with  the  thermal  instability  theory 
(sections  3.3  and  3.4).  Thus  to  reduce  instabilities,  a  high  power 
source  voltage  Uq  with  a  large  load  resistor  is  required.  From  the 
economic  point  of  view  this  is  not  ideal  since  a  large  portion  of  the 
energy  is  lost  in  the  load  resistor  R^. 

The  enhancement  of  instabilities  by  low  currents  can  also  be 
explained  by  the  thermal  instability  theory  (section  3.4).  The  exist¬ 
ence  of  other  kinds  of  column  instabilities,  however,  cannot  be  ex¬ 
cluded  on  the  basis  of  the  experimental  measurements  made.  A  more 
detailed  theoretical  description  as  well  as  direct  photographic  and 
spectroscopic  records  are  necessary  for  more  accurate  conclusions. 
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The  electrode  shape  and  material  have  a  strong  influence  on 
the  stability  of  the  discharge.  However,  the  influence  of  the  elect¬ 
rode  distance  is  slight.  On  the  basis  of  the  above  phenomena  it  is 
possible  to  deduce  that  the  observed  instabilities  are  combinations 
of  surface  and  column  instabilities.  The  surface  instabilities  occur 
in  a  short  part  of  the  plasma  column  (shorter  than  25  p) ,  but  the 
voltage  drop  in  this  space  (19  volts)  represents  an  electric  field  which 
is  higher  than  7.5  kv/cm  (previous  measurements  made  by  the  author  in¬ 
dicate  a  field  of  about  30  kv/cm) .  This  would  imply  that  even  a  small 
change  in  the  discharge  structure  in  this  space  would  result  in  a  con¬ 
siderable  change  in  the  discharge  current.  The  enhancment  of  instabilities 
at  low  currents  is  probably  due  to  a  column  instability  (thermal  in¬ 
stability)  and  is  insensitive  to  the  length  of  the  discharge  column 
because  of  the  weak  electric  field  present  in  the  column. 

The  influence  of  the  cathode  shape  on  the  discharge  instabilities 
may  be  explained  through  the  cathodic  spot  motion.  If  the  cathodic 
spot  moves  on  a  sharp  tip,  the  length  of  the  current  streamlines  in 
the  strong  electric  field  are  changed  much  more  than  if  the  spot  moved 
on  a  flat  electrode  surface.  The  electrode  material  is  probably  respon¬ 
sible  for  the  changes  in  the  electric  field  structure  in  the  proximity 
of  the  electrode  surface,  as  well  as  influencing  the  motion  of  the 
cathodic  spot. 

The  influence  of  the  liquid  used  on  the  discharge  instabilities 
is  not  clear.  The  calculation  of  the  properties  of  substances  such 
as  silicon  oil,  benzene  and  hexane  at  high  temperatures  is  a  difficult 
problem.  The  structure  of  the  plasma  column  and  the  gas  shell  as  well 
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as  the  chemical  reactions  with  the  electrode  surface  can  influence  the 
stability  of  the  discharge.  The  presence  of  carbon  atoms  in  all  three 
liquids  (hexane,  benzene,  silicon  oil)  with  small  instabilities  (see 
Fig.  4.7)  is  remarkable.  In  addition,  the  creation  of  carbides  on  the 
electrode  surface  in  EDM  in  hydrocarbons  has  been  observed.  Thus  some 
influence  by  the  carbided  electrode  surface  is  possible.  In  addition, 
the  very  strong  bonds  in  silicon  oil  cause  the  "bump"  in  the  thermal 
conductivity  (caused  by  the  diffusion  of  dissociation  energy)  to  appear 
much  closer  to  the  plasma  column  than  in  the  case  of  water.  Thus  the 
shell  of  the  hot  neutral  gas  is  smaller  than  that  in  water  and  the 
plasma  column  is  situated  in  a  "cold  tube"  where  the  possibility  for 
its  motion  is  now  more  restricted. 

5 .  3  Summary 

The  problems  which  have  been  considered  in  this  thesis  are 
closely  connected  with  practical  applications  of  plasma  discharges 
in  industry.  A  complete  and  detailed  solution  to  the  problems  has 
not  been  obtained  because  of  the  complexity  of  the  mutually  inter¬ 
acting  phenomena. 

Hov;ever,  valuable  information,  both  quantitative  and  qualita¬ 
tive  about  the  discharge  dynamics  and  stability  in  liquids  has  been 
obtained,  and  an  outline  of  the  behaviour  of  independant  generator 
discharges  has  been  given.  Techniques  for  the  possible  reduction  of 
instabilities  in  EDM  have  been  discussed.  Finally,  this  thesis  serves 
as  a  basis  for  the  further  theoretical  and  experimental  research  on 
discharges  in  liquid  dielectrics. 
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APPENDIX  A 

DESCRIPTION  OF  THE  IMPULSE  GENERATOR 

The  schematic  of  the  impulse  generator  is  given  in  Figure 
(A.l).  The  generator  is  able  to  produce  single  power  impulses  with 
approximately  rectangularly  shaped  current  impulses  for  discharges 
in  the  spark  gap  G.  The  power  circuit  consisting  of  a  power  cap¬ 
acitor  C^,  loading  resistor  R^,  switching  SCR  S^ ,  extinguishing 
capacitor  C^  and  extinguishing  SCR's  S^  and  Sg,  is  controlled  by  a 
time  delay  pulse  generator  with  SCR's  S^  and  Sg  manually  controlled 
by  the  push  botton  switch  PBS.  In  the  rest  state,  switch  PBS  is 
disconnected  and  SCR's  S^  and  Sg  are  opened.  If  PBS  is  closed, 

SCR  S^  closes,  and  capacitor  C^  is  quickly  discharged  through  the 
transformer  T  into  the  larger  capacitor  C  The  secondary  coil 

_L  1U  •  •  .  -L  -3 

of  the  transformer  T^  induces  a  sharp  impulse  which  is  damped  by 

resistor  R., _  and  used  for  closing  SCR  S..  .  Capacitor  CQ  is  simultaneously 
Id  1  o 

discharged  through  the  capacitor  C  .  Thus  the  voltage  on  the  capacitor 

O 

C  rises  according  to  the  time  constant  of  the  combination  of 

C,n  CG,  R00.  When  this  voltage  (limited  by  the  zenner  diode  D  ) 

is  high  enough,  SCR  Sg  closes  and  capacitor  C^  discharges5 creating  two 
time  delayed  impulses  in  the  secondary  coils  of  the  transformer  T^. 

After  capacitors  C-,,  Cg  and  C^  discharge  through  the  SCP.’s  S^  and  Sg , 
the  capacitors  C^,  Cg,  Cg  slowly  charge  again  through  high  resistors 
R  ? ,  R^g ,  Capacitor  Cgg  discharges  through  resistor  R^^  and  the 

variable  time  delay  circuit  returns  the  initial  state.  Time  delays 
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Fig.  (A.l)  Schematic  Diagram  for  the  impulse  generator. 
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in  the  range  of  10  psec  to  1  msec  may  be  obtained  by  switching  between 
capacitors  C-^q  ^3*  ^  regulation  of  the  time  delay  is  provided 

by  potentiometer  R  . 

In  the  power  circuit  all  SCR's  are  open  in  the  initial  state. 

The  power  capacitor  and  the  extinguishing  capacitor  C are  charged. 

The  starting  impulse  from  transformer  T1  closes  the  SCR  and  capacitor 

C  slowly  discharges  (because  large  value  of  R  )  through  R  ,  S  R 
x  j  4  15 

R6  until  the  breakdown  in  the  spark  gap  G  occurs.  The  time  de¬ 
layed  impulses  from  transformer  T  close  the  extinguishing  SCR's  S 

^  z 

and  S  and  the  extinguishing  capacitor  discharges  through  S„  ,  S„,  R 

J  2  3  11 

^2’  ^1  aS  as  as  dosed.  When  no  current  is  streaming  through 

SCR's  S^  and  S^,  these  open  spontaneously.  After  the  discharge 
extinguishes,  capacitors  C^  and  C^  charge  slowly  through  resistors  R^ 
and  R^  respectively  and  all  circuits  return  to  the  initial  state.  The 
power  SCR's  are  protected  against  overvoltage  by  avalanche  diodes 
and  .  The  non-equal  reverse  current  of  SCR's  S^  and  S^  is  compensated 
for  by  resistors  R^  and  R  .  The  difference  between  their  dynamic 
characteristics  is  compensated  by  C^,  R^.  C^ ,  R  .  The  extinguishing 
current  peak  is  limited  by  R^.  The  electric  field  which  is  necessary 
for  the  creation  of  a  graphite  bridge  is  provided  by  the  dc  voltage 
source  connected  across  the  gap  G  through  diode  D^.  Two  SCR's  are 
used  in  the  extinguishing  circuit  because  of  the  higher  inverse 
voltage  peaks  that  encountered.  Only  one  SCR  is  necessary  in  the 


main  circuit. 
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